Flow equation approach to the pairing problems 
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We apply the flow equation method for studying the fermion systems where pairing interactions 
can either trigger the BCS instability with the symmetry breaking manifested by the off-diagonal 
order parameter or lead to the gaped single particle spectrum without any symmetry breaking. We 
construct the continuous Bogoliubov transformation in a scheme resembling the renormalization 
group procedure. We further extend this continuous transformation to a case where fermion pairs 
interact with the boson field. Due to temporal quantum fluctuations the single particle excitation 
spectrum develops a gap which is centered around the renormalized boson energy. When bosons 
undergo the Bose Einstein condensation this structure evolves into the BCS spectrum. 



The flow equation method proposed a couple of years 
ago by Wegner P| and independently by Wilson and 
Glazek turned out to be a very useful theoretical 
tool in the studies of condensed matter and high energy 
physics (for a recent survey of applications see for in- 
stance |3j)- The main idea of this technique is to trans- 
form the Hamiltonian to a diagonal or to a block-diagonal 
structure using the continuous unitary transformation 
with a flexibly adjusted generating operator. 

A continuous diagonalization can be thought of as a 
process of renormalizing the coupling constants in a way 
similar to the Renormalization Group (RG) approach |j| . 
Difference between these methods is such that using the 
flow equation procedure one does not integrate out the 
high energy excitations (fast modes). Instead of this, 
they are renormalized in the initial part of the transfor- 
mation. On the other hand the low energy excitations 
(slow modes) are transformed mainly at the very end. 
In the condensed matter these low energy excitations are 
most relevant for the physical properties. 

In a first part of this paper we show how to construct 
the continuous unitary transformation for the reduced 
BCS model. Exact solution can be obtained of course 
by various ways, in particular via the Bogoliubov- Valatin 
transformation |5j] . We rederive here the rigorous solution 
using the flow equation method. This worthwhile to do 
because the pairing instability problem is in the RG stud- 
ies not a trivial issue. Only very recently there appeared 
in the literature several concepts based on the functional 
RG 0,0 discussing how to work out the Cooper instabil- 
ity. Since the flow equation method is relative to the RG 
scheme we think it would be instructive to follow Weg- 
ner et al and show how one can deal with the pairing 
interactions. 

In the second part we generalize this transformation to 
a case of fermions interacting with the boson field. Such 
situation can refer for example to the conduction band 
electrons coupled through the Andreev type scattering to 
the localized electron pairs ^ . On more general grounds 
it can also describe some boson type modes interacting 
with the correlated electrons The other realization 
of this scenario is possible in the ultracold gases where 
fermion atoms interact with the weakly bound molecules 
giving rise to the Feshbach resonance and ultimately lead- 



ing to the resonant superfluidity 



I. THE BCS PAIRING 

In order to have a simple example how the fast and 
slow energy modes are scaled during the continuous uni- 
tary transformation we consider first the exactly solvable 
bilinear Hamiltonian 



H 



h.c. 



(1) 



describing fermions coupled to some pairing field At. In 
we use the standard notation for the creation (an- 
nihilation) operator cj^^ (cko-) and energy Ck = Ck ^ M 
is measured from the chemical potential /i. Pairing field 
can be thought of, for instance, as a mean field approx- 
imation to the weak pairing potential between electrons 
Ak = - Z]q^k,q(c_qxCq|), whcrc Vk,q < 0. In general, 
one can treat Ak as the boson operator e.g. arising from 
the Hubbard Stratonovich transformation for the inter- 
acting fermion system. It can also correspond to some 
boson field responsible for mediating the pairing between 
electrons. Coupling to the boson operator will be dis- 
cussed in the next section. 

Hamiltonians of the fermion and/or boson systems 
with a bilinear structure have been already consid- 
ered using the fiow equation method P, 0, 0, E3 ■ How- 
ever, the authors have focused so far on determination 
of the quasiparticles energies. Here we supplement their 
analysis by calculating the correlation functions which 
contain information about the order parameter as well 
as a characteristics of the excitation spectrum. 

It is well known, that the bilinear Hamiltonian can 
be diagonalized by the (single step) Bogoliubov transfor- 
mation |5j. All physical quantities, static and dynamical, 
can hence be determined exactly. In this work we get the 
same rigorous solution in a process of continuous diago- 
nalization of the coupled states |k, '\> and | — k, J,>. How 
fast this can be achieved depends on a distance from the 
Fermi surface |k — ki?| (or on the energy ^k)- Simulta- 
neously with diagonalization there emerge the coherence 
factors Uk, Wk- Their final values depend on the relative 
distance of the momentum k from the Fermi surface. 
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We construct the continuous unitary transformation 
where I stands for a formal flow parameter varying 
between the initial I = value and value I > such, that 
transformed Hamiltonian H{1) — U{l)HU^^{l) simplifies 
to the needed structure (diagonal, block-diagonal, tridi- 
agonal or any other). With an increase of the flow param- 
eter I the Hamiltonian evolves according to the general 
flow equation ,,1] 



dH{l) 



dl 



m).m)] 



(2) 



where f]{l) 



dU(l] 
dl 



U ^{l). We now require the trans- 



formed Hamiltonian H(l) to preserve its bilinear struc- 
ture (Q) and let only the coupling constants to become 
Z-dependent (renormalized) 



H{1) 



const{l) 



h.c. 



(3) 



with the initial conditions ^k(O) = Ck, Ak(0) = Ak and 
const{0) = 0. This strategy reminds the scheme of RG 
approach. It can be shown Q that, roughly speaking, 
the flow parameter I corresponds to A~^. 

There are several ways for choosing the generating 
operator fi{l). In order to drive the Hamiltonian to a 
diagonal structure we adopt the Wegner's proposal [J 
m) = [Ho{l),H{l)l where HoU) = Ek,.a(OCck. de- 
notes a diagonal part. With this choice the Hamiltonian 
becomes diagonal in the limit Z — > cxo. In explicit form 
the generating operator for J^Jl reads 



k 

Substituted into the equation Q it gives 



(4) 



dl 



'^''^^^ ^ 4^a(0iAk(0P(4.ek.-i) 

k,cr 

- 4^(a(0)'(Ak(o4TeUi + /^-' 



which is identical with the following set of flow equations 



dl 
dl 



4^k(0|Ak(0h 

-4(ek(0)'A;;(/) 



(6) 
(7) 



and dconst{l)/dl = -4^k,o-Ck(0|Ak(0P- 

Formally, solution of the equation ||7J) can be given as 



|Ak(/)| 



lAwIe" 



'4 J^dl'[(^il')f 



(8) 



which yields that lim;^ooAk(/) = for all k 7^ kp. 
When we multiply equation © by ^k(0 ^^^^ equation {Tj) 
by Ak(0* then their sum gives the following invariance 



dl 



{(a(0)' + |Ak(OP}=0. 



(9) 
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FIG. 1: Renormalization of Ak(0 and (,k{l) during the flow 
for several values of I as denoted in the legend. At the start- 
ing point Z = we used the isotropic coupling Ak(0) = A. 
Energies are expressed in units of A while I in units of A~^. 



This constraint together with vanishing of Ak(/) in the 
limit I —>■ 00 implies the known Bogoliubov spectrum 



^k(cx)) = sgn(^k)A/(|k 



(10) 



In figure ^ we plot the pairing field Ak(0 and the 
change (renormalization) of fermion energies ^k(0 ~ Ck 
at several stages of the continuous transformation. We 
notice, that fast modes (states distant from the Fermi 
energy) are transformed rather in the first part of the 
process and change of their energies is rather small. The 
slow modes (i.e. the low energy excitations) have to be 
worked on much longer. Asymptotically, at Z — > 00, the 
entire spectrum reduces to the Bogoliubov structure i{TU\i . 
Now we turn attention to the dynamical quantities 



(5) which can be expressed via the normal ((ck 



and the anomalous single particle Green's function 
((ck|; c_kx))a)- We introduced the standard Fourier 
transforms for the retarded fermion Green's function 
/ duje"^mA;B))^ = -~iQ{t){A{t)B + BA{t)) and time 

evolution is given by A{t) = e**^Ae~**^. 

Thermal averaging for some arbitrary observable O 
is defined as (O) = Tr je-'^^oj /Tr je-^*^!, where 

/3~^ ~ ksT. Since the trace is invariant under unitary 
transformations we can write down 



Tr 



Tr 



Tr 



i^Vl{l)e-^^dU-\l)} 
{e-^^«6(0}, 



(11) 



where 0(l) — U(l)OU^^{l). It is convenient to com- 
pute such trace in the limit I = 00 because Hamilto- 
nian becomes then diagonal. However, the price which 



3 



we pay for having the simplified Hamiltonian at / = oo 
is the additional necessity to transform the observables 
O 0{l) 0{oo). This is to be done similarly as with 
the Hamiltonian Q using 



dOil) 
dl 



[77(0,0(0]- 



(12) 



With the choice of the generating operator 77(0 given 
by Q we find from the commutator [7)(/),Ckf]) that for 
I > the annihilation operator Ckf gets convoluted with 
cLkj- It is natural to propose the following Ansatz for 
the /-dependent operator 



CkT(0 = '«k(OckT + t'k(OcLkj,- 
Similar analysis for c^^.^ operator leads to 

c-ki(0^ = -'fk(OckT +Wk(Oclki 



(13) 



(14) 



with the initial conditions Uk(0) — 1, Wk(0) = 0. These 
parametrized equations H13I14|I substituted into 1)12(1 with 
r]{l) given in Q lead to the coupled flow equations 



dl 
dl 



2a(OAk(0'^k(o, 

-2ek(0Ak(0^^k(/). 



(15) 
(16) 



After straightforward algebra using equations H15llt)|) 
we obtain the following invariance |7ik(0P + |''^k(OP — 
const = 1 which assures that the /-dependent fermion 
operators (|13I14|I preserve the anticommutation relations 



|ck<T(07Ck'<T'(0} = (^k.k'^CT 



(17) 



Without a loss of generality we assume Ak(0 to be real 
so that the unknown coefficients Mk(0, ■'^k(0 become real 
too. To determine their asymptotic I = 00 values we can 
rewrite ()15|l as 



du]^il) 
«k(0 



2ek(0Ak(0'^^ 



(18) 



/— 00 



and we further integrate both sides in the limits /;_g 
Using i;k(0 — ^1 — (uk(0)^ S^t for the l.h.s. 

■ °° rfuk(0 



1=0 - iMi)y 



arcos Uk(oo) 



(19) 



due to Uk(0) = 1. Using equation iQ we can replace 
2^k(0Ak(0'^^ by ~'^Ak(/)/2^k(0 from the invariance 
© we obtain that ^k(0 = V^li^) - |Ak(OP- The r.h.s. 
of equation ((TH|l gives after integration 
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FIG. 2: Flow of the coherence factor ul.{l) for three fermion 
states of the initial energy = 0.01, 0.1 and 1 (in units A). 
Inset shows the effective value of the coherence factor u^(oo) 
versus ^k. 



because Ak(oo) — 0. Combining the results H19I2U|) 



f ^ |Ak| 
2 arcos Wk(oo) = -- — arcos j——, — r- 
2 l?k(oo) 



we finally determine the I — 00 factors 



"k(oo) 



Ck(oo) 



i-«kM- 



(21) 



(22) 



One can simply show that 2Mk(oo)wk(oo) — Ak/|^k(oo)|- 
In figure El we present the coherence factor Wk(0 a 
function of the flow parameter /. Again we notice that 
for momenta distant from the Fermi surface the coher- 
ence factors evolve rather quickly from the initial value 
M^(0) = 1 to the asymptotic values (see the inset). The 
coherence factors of the slow modes establish later on, 
similarly to renormalization of the corresponding excita- 
tion energies shown in figure^ Saturation occurs around 

I cx i/a(oo)'- 

Since the transformed Hamiltonian iJ(oo) is diagonal 
we can easily derive the single particle Green's functions 
(and the higher order Green's functions too). With the 
parameterizations (|13I14|I we obtain 



((ckT;ci,|)). 
((ckT;c-ki)). 



■ ?k(oo) 



.... , a;-|-^k(oo) 
Mk(oo)t^k(oo) _ ^^k(oo)^^k(oo) 



- ?k(oo) 



UJ - 



-Sk(oo) 



(23) 
.(24) 



The equal time expectation values defined by (BA) = 
-i/:°^dL^/(c^,r)Imag((i;S))^+,o+ (where /(l^,T) = 
^pPuj _|_ 2j ^ is the Fermi distribution function) yield the 



equations for average momentum occupancy 



^ktCkT) = 1^ 



_ 4L_tanh^' 
^k(^) 2kBT 



(25) 
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and for the ofF-diagonal order parameter 



(C-kiCkt) = 



tanh 



Ck(oo) 



2ek(w) 2kBT 



(26) 



These equations (|25I26|) exactly reproduce the rigorous 
solution of the reduced BCS Hamiltonian. 

In appendix A we explain how to generalize the present 
treatment to account for the scattering of finite momen- 
tum fermion pairs. 



II. COUPLING TO THE BOSON MODE 

We apply here the same method to the non-trivial 
problem describing itinerant fermions coupled to some 
dispcrsionless boson mode 



H 



k.o- 



kCT 



Cko 



9k 



kT^'-ki 



h.c. 



(27) 



Here the boson mode can be regarded for instance as 
a pairing field derived from the Hubbard Stratonovich 
transformation in the system of interacting fermions. 
Model (|27|l can describe also the Andreev tunneling be- 
tween c - fermions and some pair reservoir denoted by 6 - 
particles. Other possibility is to think of the case where 
itinerant fermions (e.g. conduction band electrons) coex- 
ist and interact with some localized fermion pairs 9] . The 
model (|27l) is also often applied to describe the ultracold 
fermion atoms coupled to the weakly bound molecules 
effectively leading to resonant Feshbach scattering 

To keep a conserved total number of particles Ntot = 
Ek (T ^ka^ko- + 26^6 we apply the grand canonical ensem- 
ble, stands for the boson energy measured from 2/i and 
(7k denotes the coupling constant. It can be shown that 
physical properties of this model depend only on a mag- 
nitude of (7k J in other words all |(?k|e*'^'' lead to identical 
results independently of phase (jj^. 

In the mean field treatment of the Hamiltonian H27|l 
one usually introduces the linearization bcl^-^c^_^^^ — 

(6)cj^|C^i^j^ + 6(cj^|'5^ki) ^9]. Such idea is based on as- 
sumption that there exists a finite amount of the Bose 
Einstein (BE) condensation of b particles. Hamiltonian 
(|27|1 simplifies then to the ordinary BCS problem, where 
A = —g{b). Let us however remark that the BE conden- 
sation of infinitely heavy boson (characterized by a dis- 
crete energy level il) cannot occur. Our present analysis 
based on the fiow equation procedure shows a possible 
route to go beyond such mean field approximation. 

We construct the continuous canonical transformation 
which decouples c from b particles. This can be achieved 
only approximately by choosing the following generating 
operator 



m^ii) - m] 9k{i)bciAki ~ ^.c. (28) 



Applying (|28|l to the flow equation ^ and going through 
appropriate normal ordering we obtain the /-dependent 
Hamiltonian 



H{1) ^a(/)4.^k. + E(5k(0^4Tcl 



k,a k 

+ fi(/)&^& + $]C/k.k'(04TcUiC-kuCk'T (29) 

k,k' 

with [/k,k'(0) = 0. In a straightforward way we derive 
the following set of flow equations 



dl 

dn{i) 
dl 

dg^il) 
dl 

dUky{l) 
dl 



= 2|(7k(on2^k(o-f^(0]- 



(30) 



= 2Y,\9k{l)? [2Ul)~m\ (2nL-l][31) 



= -gk{l)[2£,^{l)-m] 

= -gk9l'{l)[^Ul)-m\ 
-gk'9Ul) i2^k'il)~m))- 



(32) 



(33) 



Here n^^ denotes the average fermion occupancy of the 
state |k, (t) and ~ (^^^}- From formal solution of the 
equation ^ 5k(0 = ffk exp{- dl'[2(i,{l') - n{l')]^} 
we notice that g]s_{l—>oo) — 0. 

In order to close the set of flow equations H3UI30|I we 
must determine the distribution function n^^. From 
analysis of the flow equation H12I) for the annihilation 
and creation fermion operators we come to a conclusion 
that 

CkT(0 = tik(0 CkT + «k(0 ^ (34) 
c^_ki(0 = -«k(0 CkT +Mk(0 cLki- (35) 

Here the Z-dependent factors satisfy the following flow 
equations 



Ukjl) 

dl 
dl 



= v^{l)gu{l) [2£k(/) - m] (n'' + n^ki) (36) 
= -uuil)gk{l)[2eu{l)-m]- (37) 



with the boundary condition Wk(0) = 1 and Uk(0) = 0. 
For simplicity in (|36I37|) we neglected the terms pro- 
portional to (7k(0t^ki,k2(0 which might eventually con- 
tribute some higher order corrections of the order ~ g^. 
Combining (|36|l and 137(1 we obtain the invariance 



|wk(OP + (n^ + "-^ki)bk(OP ^ const 1 



(38) 



which guarantees that operators c^^W defined in (j^ 
obey the anticommutation relations. 

Various expectation values (O) are easy to carry out in 
the limit I — > oo because fermions are decoupled from the 
boson field. Some delicate problem causes the two-body 
interaction C/k,k'(oo)- As will be shown below this poten- 
tial is small, so in the lowest order perturbation theory we 
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FIG. 3: The single particle excitation energies (top panel) 
and their corresponding spectral weights (bottom panel) . The 
thin dashed line indicates that the excitations Bk.i/ are gaped 
around the renormalized boson energy if2(oo). We used for 
computations 5k = 0.05-D, 0/2 + /i = -0.3D, ksT = 0.015-D. 



incorporate its effect via the Hartree shift ?7k,k(oo)ri^^ to 
fermion energies. This weak point of the present analysis 
can be however systematically improved. 

Using the parameterization (|34I35I) we find the normal 
single particle Green's function 



((ckT;c[T))'' 



|uk(oo)P 



+ (n^ + ; 



LU - i?k.2 



(39) 



with quasiparticle energies E'k.i = Ck(oo) + UuM{oo)n^^ 
and £'k.2 — ^^(oo) — E'k.i- The corresponding spectral 
function A{'k,uj) = — ilmag((ck|; cjj.|))(j+io satisfies the 
sum rule /^(iwA(k, = 1 due to invariance (|38|l . Fi- 
nally, the momentum distribution function is thus given 

by 



|Mk(oo)P 



+ (n^+^ 



|uk(oo) 



1 



(40) 



We solved numerically the set of coupled flow equa- 
tions H30I33I36I37|) and ^ using the Runge Kutta al- 
gorithm. In our calculations we considered the finite 
fermion band —D/2 < < D/2 and assumed the flat 
density of states. Total particle concentration was fixed 
for (Ntot) = 1 which corresponds to nearly equal popu- 
lations of fermions and bosons n^^ ^ . 

In the top panel of figure Q we show the single parti- 
cle fermion spectrum which turns out to be gaped. Two 
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FIG. 4: Potential f/k,k(oo) of the two-body interaction be- 
tween fermions for the same set of parameters as in figure 
3. 



branches of the excitations £'k,i/ are discontinuous around 
the energy ri(cxD)/2 (the dashed line). For here consid- 
ered situation it is located slightly above the chemical po- 
tential 11 and is dependent on temperature. We checked 
that the value of the gap is rather independent of T . 
The corresponding spectral weights are shown in the bot- 
tom panel of figure 13 They behave similar to the BCS 
formfactors but we notice only a negligible dependence 
on temperature. Figure 01 additionally confirms that the 
two-body potential is indeed weak. 

We would like to emphasize, that the gaped spectrum 
shown in figure is not related to any off-diagonal order 
parameter. According to the parameterization H34I I35|) 
we easily determine that the anomalous Green's function 
((clkj,! '-kt)) identically zero. The above mentioned 
structure should hence be referred as a pseudogap. The 
broken symmetry can arise if and only if a certain fraction 
of h particles becomes BE condensed . There are two 
possibilities for this to occur: a) either we assume bosons 
to be mobile right from the outset of the problem, or b) 
we allow for a finite momentum exchange q ^ in the 
interaction term c\^^':^„_\^\b<\- 

The second option has been explored in the literature 
by several groups using various many-body techniques 
(for representative list of the references see the review 
paper jlTj). The purpose of our present study was to 
show that the temporal quantum fluctuations alone can 
induce the (pseudo)gap structure - there is no need for 
involving spatial fluctuations. Of course in real systems 
there always exist both, temporal and spatial fluctua- 
tions. The latter cause that upon lowering temperature 
the (pseudo)gap smoothly evolves into the true gap of 
superconducting state (l8| . 



III. CONCLUSIONS 

By means of the flow equation method 0] we ana- 
lyzed the bilinear Hamiltonians describing fermion sys- 
tems with pairing interactions. This new technique be- 
comes more and more popular Q owing to its conceptual 
simplicity which allows to go beyond the frame of various 
perturbative approximations |l9l |. This method is based 
on the continuous canonical transformations in course of 
which all the parameters are gradually renormalized. In 
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some way this can be compared to renormalization of 
the couphng constants within the RG procedure 3 . The 
high energy excitations are here renormahzed mainly in 
the beginning, while the low energy excitations must be 
worked on until the end of transformation. 

First part of our paper illustrates how to reproduce 
the rigorous solution of the BCS model using a contin- 
uous version of the Bogoliubov transformation In 
particular, we show how one can approach the symme- 
try broken problems which usually is not an easy task 
for the RG methods |^ 0] ■ This scheme can be extended 
for analysis of the two-body interactions where in princi- 
ple various kinds of instabilities can arise. Some results 
for the 2-dimensional Hubbard model have been already 
reported with use of the flow equation method. 

In the second part we focused on the case, where 
fermion pairs are coupled to some single level boson field. 
Such situation can take place when the correlated fermion 
system is afi'ected by bosonic modes such as e.g. the pair- 
ing fluctuations in the systems of reduced dimensionality 
(for instance the HTS cuprates) Moreover, this sort 
of physics is recently intensively studied for the ultra- 
cold fermion atoms where interaction with the weakly 
bound molecules leads to the Feshbach resonance which 
is a driving mechanism for the atomic superfluidity p3 |. 

We solved selfconsistently the corresponding set of flow 
equations with a help of the numerics. We showed that 
interactions are responsible for appearance of a gap in the 
single particle fermion spectrum. This gap is centered 
around the boson energy ^n{(x) instead of the chemi- 
cal potential. It does not signify any symmetry break- 
ing because no order parameter is present in the system. 
The off-diagonal order parameter can eventually appear 
if bosons have a finite mobility (characterized by some 
dispersion f2q different from the discrete energy n con- 
sidered here) and if those bosons undergo the BE conden- 
sation jl^ . The center of gap moves then to the chemi- 
cal potential as has been explained in the previous study 
|l5j . Here, we emphasize that the single particle gap 
can appear in a normal state purely due to the tempo- 
ral quantum fiuctuations. In realistic systems with addi- 
tional spatial fluctuations such normal state (pseudo)gap 
is expected to evolve smoothly into the gap of symmetry 
broken superconducting state. 
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APPENDIX A 



Hamiltonian 



k,q 



4T<-ki + /i-c-) (Al) 



which reduces to ^ if Gi^q — (5q_oAk. Besides q=0 this 
model allows also for a scattering of the finite momentum 
q fermion pairs. In what follows below we briefiy analyze 
the effect of such finite momentum scattering. 
The generator from Wegner's proposal 



V 



'(0 = -EKk(0 + eq-k(0](Gk,q(04T<-ki)-^-C- 



k,q 



(A2) 



applied to the flow equation ((2Jl induces the off-diagonal 
terms cj^^Cp^ko-- In order to eliminate them from the 
transformed Hamiltonian H{1) we replace (|A2|) by 



m = v"ii) + E ^k,q,.(o (eLsq- - C^.Ck.) • (A3) 

k,q,a" 

After simple algebraic calculations we flnd that the off- 
diagonal terms are exactly canceled if 



7k.q,T(0 ^ 



^k(0 + ^k'(0 



G'k',k'+q-k(OGk,q(Oj 



(A4) 



7k 



.qa(0 = E 



ek(o+ek'(o 



Gk',k+k' (OGk',k'+q(0- 



(A5) 



The modified generating operator (|A3fl has a virtue to 
preserve the initial structure of the Hamiltonian H{1) = 

Ek,^^k(OcLck^ - Ek,q (G'k,q(04TCq-ki + 'i-C-) + 

const{l). The corresponding set of fiow equations is 



dl 



= 2E(a(0+eq-k(0)|Gk,q(0l', (A6) 



'^^^ = - m) + eq-k(O)' G^,q(0 (A7) 
-I-2E (7k,k',<T(0 - 7k',k,<T(0) Gk',q-k+k'(0> 



d const{l) 
dl 



-2E(^k(0+eq-k(0)|Gk,q(0l' 
k,q,fT 



(A8) 



We were not able to solve analytically the equations 
(|A6IA8|l therefore we explored them numerically assum- 
ing the tight binding dispersion = — 2t cos{kxa) and 
using the pairing potential Gk.q in a Lorentian form 



The continuous Bogoliubov transformation (discussed 
in section 1) can be extended to the more general bilinear 



Gk,q = AA(n) 



A 



n [q.^a] + 1 



(A9) 



The normalization factor M{n) was taken such that 
EqG'k,q = A and we set A/4i = 0.01. 

Flow of the operators c^^il) takes here the form 

CkT(0 Xl"k,q(OCq+kT +Xl^k,q(OCq_ki (AlO) 



(0 



XI "k,q(OCq+kT + XI "k,q(OCq_l(f 11) 



which generalizes the previous equations I14|l due 
to finite momentum scattering. Substituting (|A10IA11|) 
into the flow equation H12|l we obtain 



(iMk,q(/) 



= Kk(0+eq-k(0]Gk,q(0'^k,o(/) 



+ EKk-qa)+eq+p(0]Gq-k ,2q-k+p 

P#k 

(A12) 

= - Kk(0 + eq-k(0] Gk,q(0^k.o(0 

- X Kk-q(0 + Cq+p(0] Gq_k,2q-k+p(0^k,q(0- 
P#k 

(A13) 



From these equations (jA12IA13|) we derive the invariance 

El"k,q(OP+ElWOP = l (A14) 




1 



CO/A 



FIG. 5: The local density of states p{ijj) — ^^Aiyi,ui) of 
fermions coupled to the potential Gk.q given in <A9ll with 
n = 10, 50 and 100. The density is normalized to Po('^) of 
the nointeracting sytem. 




which again leads to the anticommutation relations (|17|l . 

With use of the parameterization (|A10IA11(I we can 
now determine the single particle Green's functions. 
Let us first consider the diagonal part (in the Nambu 
representation). The spectral function j4(k, cj) = 

— 7r^^Imag((ck| ; cj^|))a;+jo+ turns out to consist of two 
contributions 



>l(k, Uj) = Acoh{K ^) + ^mc(k, Uj). 

The coherent part 



(A15) 



Acoh(k,w) = |Mk,o(oo)p (5[cj -,^k(oo)] 

+ |i;k,o(oo)p ,5[L^ + e-k(oo)] (A16) 

describes the long-lived quasiparticle modes of energies 
i^k(oo) similar to the Bogoliubov modes pU|) discussed 
previously for Gk.q = Ak(5q,o. The remaining incoherent 
part 

AinciK^) = |Mk,q(cx3)P 5 [u - ^q-|-k(oo)] 

q#0 

+ X l^k,q(oo)|2 ,5[a; + eq-k(oo)](Al7) 

q#0 

corresponds to the background spectrum which is spread 
over the large energy regime. Due to invariance (|A14|) 



FIG. 6: The single particle spectral function yl(k,tj) obtained 
for the potential 1A9II with n — 100. The thick vertical lines 
represent the spectral weights |uk,o(oo)P and |uk,o(oo)p of 
the long-lived quasiparticles ±^k(oo). The thin solid lines 
show UJ dependence of the incoherent background yli„c(k,aj). 



the total spectral weight of the coherent and incoherent 
parts satisfies the sum rule dujA^k^uo) — 1. 

Figure 13 presents the local density of states calculated 
by us with use of the definition p{u}) = A(k, w). We 
notice that the excitation spectrum is not longer truly 
gaped. There occurs only a partial suppression of the 
fermion states around the energy [u — Q). This property 
should be assigned to the scattering of finite momentum 
fermion pairs on the potential Gk,q introduced ad hoc in 
(|A1|) . Physically this means that only a certain fraction 
of the fermion states is expelled from a vicinity of the 
Fermi surface. 

Figure lO shows both parts of the spectral function 
for the pairing potential l|A9|l with n= 100, which yields 
the strongest suppression of the low lying fermion states 
ploted in the figure [S] Two modes of the coherent 
part have similar spectral weights as the BCS coher- 
ence factors. However, in distinction from the standard 
BCS solution, the corresponding quasiparticle dispersion 
±^k(oo) is in this case gapless. This is evidently caused 
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by a coupling to the finite momentum fermion pairs. As 
regards the incoherent background it builds up mainly 
near the Fermi surface hp- Such damped fermion states 
are located around the quasiparticle energies ±fk(oo)- 
We estimated that for k = k^? the incoherent states con- 
tribute nearly 30 percent of the total spectral weight. 

The off-diagonal single particle Green's function 
((ck|; C-pi))u: has a similar structure to IjAlTp . Skipping 
the unnecessary technical details we show the expression 
for expectation value of the order parameter 



(C-piCkt) 



1 \p f Mk,q+p-k(oo) •i^p,q(oo) 

2^1 exp[/3eq+pM] + l 



Mp,q(oo) t^k,q-p+k(oo) 

exp [-/3^q„p(oo)] + 1 



(A18) 



We investigated numerically the q-dependence of the or- 
der parameter (cq_kj.Ck;f) and we noticed its rather fast 
decrease against |q|. For n — 10, 50 and 100 we obtained 
that magnitude of the order parameter at |q| = 0.27r/a 
is circa 100 times smaller than (c-kxCkt)- This indicates 
that finite momentum fermion pairs are less favored in 
the system. Moreover, the finite momentum fermion 
pairs are damped entities because they arise from the 
incoherent part of the off-diagonal spectral function ex- 
isting in a broad energy interval lo for momenta ~ kj?. 
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